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Abstract 

When we have two expansions of physical quantity around two different points in 
parameter space, we can usually construct a family of functions, which interpolates the 
both expansions. In this paper we study analytic structures of such interpolating func¬ 
tions and discuss their physical implications. We propose that the analytic structures 
of the interpolating functions provide information on analytic property and Stokes phe¬ 
nomena of the physical quantity, which we approximate by the interpolating functions. 
We explicitly check our proposal for partition functions of zero-dimensional theory 
and Sine-Gordon model. In the zero dimensional Sine-Gordon model, we compare our 
result with a recent result from resurgence analysis. We also comment on construction 
of interpolating function in Borel plane. 
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1 Introduction 

In most non-exactly solvable problems, one tries to look for a small parameter or a large 
parameter so that we can set up a perturbative expansion. This perturbation series is typi¬ 
cally an asymptotic series which is divergent and has zero radius of convergence. Often the 
perturbation series takes different forms depending on the argument of its expansion param¬ 
eter. This behavior is known as the Stokes phenomenon [1]. For example, if we consider 
a physical quantity having an integral representation as in quantum held theory and string 
theory, then the integrand usually has multiple saddle points and its perturbative expansion 
with the parameter g typically takes the form 

k k k 
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where Si{g) is the “action” evaluated at each saddle point and each sum denotes the small- 
expansion around that saddle point. In standard situations, the weight is exponentially 

suppressed for real positive g and vanishes as approaching g -|-0. In other words all the 
terms except the first one in (1.1) describe non-perturbative effects. However, if we change 
8iTg{g), then Si{g) might have negative real part. For this case, the weight is not 

exponentially suppressed but exponentially growing. This type of transition happens when 
we cross an anti-Stokes line given by Re(S'j) = 0. On the other hand, the coefficients 
themselves might jump as we change arg( 5 f). This happens when we cross a Stokes line defined 
by Im(S'j) = Im(S'j). Thus Stokes phenomena has deep connections to non-perturbative 
effects. 

In most practical situations, we can access to only first few terms of a perturbation series 
around single saddle point for particular argument of expansion parameter. For this case, 
it is hard to find when Stokes phenomena occur and when contributions from other saddle 
points become important. In this paper we develop a tool to study Stokes phenomena in 
somewhat special situations. 

Sometimes we can have two perturbative expansions of physical quantity around two dif¬ 
ferent points in parameter space, for example, in theory with S-duality, field theory with 
gravity dual, lattice gauge theory with weak and strong coupling expansions, statistical sys¬ 
tem with high and low temperature expansions, and so on. One can then construct functions, 
which interpolate these two expansions. The most standard approach to this is (two-point) 
Fade approximant, which is a rational function having the two expansions up to some orders. 
Recently Sen has considered another type of interpolating function, which has the form of a 
Fractional Power of Polynomial (FPP) [2]. After a while, one of us has constructed a more 
general class of interpolating functions described by Fractional Powers of Rational function 
(FPR) [3], which includes the Pade approximant and FPP as special cases. It has turned 
out that these interpolating functions usually provide better approximations than each per¬ 
turbative expansion in intermediate regime of the parameter, see [2, 3, 4, 5, 6, 7] for various 
applications^. Although these are quite nice as first attempts, properties of interpolating 
functions themselves have not been extensively studied yet. 

In this paper we propose new properties of the interpolating functions. We focus on 
analytic structure of the interpolating functions^ by treating them as complex functions and 
discuss their physical implications. The FPR is some power of a rational function and the 
rational function has poles and zeros. When the power is not integer, then the poles and zeros 
of the rational function give rise to branch cuts of the FPR. Here we propose that the branch 
cuts of the FPR encode information about the analytic property and Stokes phenomena of 
the physical quantity, which we try to approximate. More concretely we propose that each 
branch cut of the FPR has the following possible interpretations. 

^ There is also another type of interpolating functions [8], which is not special case of the FPR. This has 
been applied to 0{N) non-linear sigma model. 

^ Note that our interpolating function approach is different from resurgence approach, which has been 
recently studied in a series of works [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]. While we use single 
types of weak and strong coupling expansions as input data, the resurgence approach uses weak coupling 
expansions around multiple saddle points and does not use strong coupling expansion. 
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1. The branch cut is particular to the FPR and the artifact of the approximation. Namely, 
this type of branch cut is not useful to extract physical information. 

2. The physical quantity, which we approximate by the FPR, has a branch cut near from 
the branch cut of the FPR. Namely, the branch cut of the FPR well approximates the 
“true” branch cut of the physical quantity. 

3. Near the branch cut, one of perturbation series of the physical quantity changes its 
dominant part. This case is further separated into the following two possibilities. 

(a) We have an anti-Stokes line of the perturbative expansion near the branch cut. 
Namely, although the perturbative series itself does not change its own form, 
contributions from other saddle points become dominant across the line. This 
possibility likely occurs for first branch cut measured from a specific axis where 
we construct interpolating functions. 

(b) The perturbative series itself does change the form. Namely, we have a Stokes line 
near the branch cut, whose diagonal multiplier is different from 1. When we have 
a Stokes fine across which sub-dominant parts of the perturbative series change, 
the FPR cannot detect this type of Stokes line. 

We explicitly check our proposal in two examples: partition functions of the theory 
and the Sine-Gordon model in zero dimensions. Similar features seem to appear also in other 
examples such as BPS Wilson loop in 4d W = 4 Super Yang-Mills theory, energy spectrum 
in Id anharmonic oscillator etc [24]. We expect that our result is applicable in more practical 
problems, where we do not know exact solutions. One possible utility of such an analysis 
is that we can anticipate analytic property and Stokes phenomena of physical quantity by 
looking at analytic structures of interpolating functions. 

This paper is organized as follows. In section 2 we introduce our interpolating functions 
described by the fractional power of rational functions (FPRs). In section 3 we study inter¬ 
polating problem in the Od theory in great detail. We look at the analytic property of 
the interpolating function as a complex function and propose its physical interpretation. In 
section 4 we analyze the Od Sine-Gordon model and check that our proposal is true also for 
this model. Section 5 is devoted to conclusion and discussions. In appendix A we compare 
our interpolating function with a recent result from resurgence analysis [20] in the Od Sine- 
Gordon model. Our result implies that the FPR and resurgence play complementary role 
with each other. In appendix B we explain an attempt to construct interpolating function in 
Borel plane and test its utility in the Od Sine-Gordon model. In appendix G we write down 
explicit forms for interpolating functions used in the main text. 


3 


2 Interpolating function 


We introduce the interpolating functions in this section, which is essentially a review of 
[3]. Suppose that we wish to determine a function F{g), which has^ the small -(7 expansion 
Fs^'‘\g) and large-^f expansion Fi^‘\g) taking the forms 

N, Ni 

= = (2.1) 

fc =0 fc =0 

We can then naively expect that these expansions approximate F(g) as 

F{g) = = F[^^\g) + Oig”-^^-^). (2.2) 


Although this seems to be a somewhat limited case, this situation includes a large class of 
physical problems, e.g., theory with S-duality, field theory with gravity dual, lattice gauge 
theory with weak and strong coupling expansions, statistical system with high and low tem¬ 
perature expansions, and etc. 

In terms of the two expansions, one can construct the following function [3] 



1 + ELi dk9^ 


(2.3) 


where 

If ^ a — b\ If ^ a — b\ 

P = - [m + n + l - , q = - [ m + n + 1 -. (2.4) 

2 \ a J 2 \ a J 

The coefficients Ck and dk are determined such that series expansions of Fm)i(g) around 5 ^ = 0 
and g = oo reproduce the small-(yf and large-^r expansions (2.1) of F{g) up to and 

0(g’^~'^~^), respectively. Due to this property, the function Fmjiig) interpolates the small- 
g and large-g- expansions up to these orders. Since the interpolating function is usually^ 
described by the Fractional Power of Rational function, we call this type of the interpolating 
function “FPR”. Note that the rational function inside the square bracket in (2.3) is a ratio 
of polynomials, he., 

p,q E Z>o, (2-5) 


which leads the following condition 


a = 


a—b 

2^+1 

a—b 

2i 


for 

for 


m + n : even 
m -|- n : odd ’ 


with £ E X. 


( 2 . 6 ) 


It is now easy to see that the FPR includes the Fade approximant and the Fractional 
Power of Polynomial (FPP) constructed in [2] as special cases. If we take 2^ + 1 = a — 6 for 
a — 6 G Z and m + n to be even, then this becomes the Fade approximant while if we take 
2^+l = m + n + l {2£ = m + n + 1) for even (odd) m + n then we get the FPP. Therefore 
we below refer to also the Fade and FPP as FPR. 

^ More generally we might have perturbative expansions around g = gi and g = g 2 with g 2 > gi and would 
like to construct their interpolating functions. However, if we change the variable as x = {g — gi)/{g 2 — g)^ 
then this problem is reduced to interpolating problem of small-x and large-x expansions. Thus our setup 
does not lose generality in this sense. 

^ When a — b is irrational number, the power a is irrational number. 
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Figure 1: Summary of Stokes phenomena for the small-^- expansion in the partition function 
of the Od theory. The blue solid lines denote the Stokes lines while the green dashed lines 
denote the anti-Stokes lines. 


3 Partition function of zero-dimensional theory 


In this section we study interpolation problem for the partition function of the Od theory. 
Although this example has been already studied well in [2, 3] for real non-negative coupling 
constant, here we consider general complex coupling. As mentioned above, we study analytic 
properties of interpolating functions and their physical implications. Let us consider the 
integral 


n9) = 


1 r 

V9 J-c 


dx e 29 ' 


(3.1) 


which can be exactly performed as 


1 

7re32fif^ 


F(g) = 


is 


1 

7re32fif^ 


/ 1 ( 


V32A 


V 32g2 


4V-g2 L 


(3232) +di 


for Re((yf) > 0 
for Re{g) < 0 


(3.2) 


where I{z) is the modified Bessel function of the first kind. The small-^r expansion of F[g) 
depends on the argument of g, and its dependence is given by 

( - ‘i\f^g^ -I- -F 0{g^) for Sirg{g) E (- 

^( 9 ) = I (105^^ -6g'^ + 2 + 0{g^)) - ^/Fie'^ (105/ -F 6/ -F 2 -F C>(/)) for avg{g) G (f, vr) 

V \f\ (105/ — + 2-1- 0{g^)) + (105/ + 6 g‘^ + 2 + 0{g^)) for arg( 5 ') G (—vr, —|) 

(3.3) 

This dependence on arg( 5 f) clearly reflects the fact that we have Stokes lines oriented along 
arg( 5 f) = ±7r/2 and tt, across which the form of the small-(yf expansion changes. Notice that 
across the ray |arg(p)| > 37r/4, the terms involving the exponential factor (in line 2 and 3 
in eq. (3.3)) are dominant compared to those without the exponential factor. This indicates 
that we have the anti-Stokes lines oriented along arg( 5 ') = ±37r/4. 

It is easy to understand this behavior from the standard saddle point analysis for \g\ <C 1. 
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Saddle points x* of the integration are given by 


X* = 0,x±, with x± = ±——. 

2v^ 

Then, the “action” S'(x) = |^ + x^ at the saddle points x* takes the values 

S'(x* = 0) = 0, S'(x* = x±) =(3.5) 

For |arg( 5 ')| < 7r/2, we can pick up only the trivial saddle point x* = 0 by deforming the 
original integral contour (—oo, oo) to a steepest descent contour, while we can pick up all the 
saddle points otherwise. We have relative minus sign in the contributions from the non-trivial 
saddle points x* = x± because directions of the steepest descent through x* = x± are opposite 
between the cases for 7r/2 < arg(p) < tt and —tt < arg(p) < —7r/2. Note that the real parts 
of the action at the non-trivial saddle points change their signs across SiTg{g) = ±37r/4. This 
means that we have anti-Stokes lines of the small-^f expansion at a-'rg{g) = ±37r/4. This 
discussion is summarized in hg. 1. The large-expansion, on the other hand, is independent^ 
of arg( 5 f): 

F(s) = (^r(i/4) + lr(-i/4)9-' + ^r(i/4)j-2 + o(j-’)) . (3.6) 



3.1 Interpolation along positive real cixis 

Let us first take the coupling g to be real and positive as usual. For this case, we have the 
following small-p and large-p expansions 


= «2‘+i = 0. (3,7) 

k=0 

= (3-8) 

1=0 • V / 

which are compared with the exact result in fig. 2 [Left]. In terms of these expansions, we 
can construct FPR-type interpolating function Fm)i{g) (see app. C.l for explicit forms). In 
fig. 2 [Right], we test validity of the FPRs by plotting 

Ftl - F(cj) 

F{9) ’ 


against g for some We easily see that these interpolating functions provide good 

approximations to the original function F{g). Especially FQ^Q‘^\g) approximates the exact 

^These behaviors of the expansions can be understood also from viewpoint of differential equation for 
the Bessel function although we can know this information after finding the exact result. The modified 
Bessel function has essential singularity at ^ = 0 and has only a branch cut singularity at ^ = oo. The 
Stokes phenomenon of the small-^ expansion is a manifestation of the essential singularity. There is no such 
behavior near the branch point, and the large g expansion is independent of arg(^). 
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Figure 2: [Left] The partition function of the Od theory (black solid), its small-^- ex¬ 
pansions Fs^^\g) (red dotted) and large-p expansions (blue dashed) for a,rg{g) = 0. 

[Right] Comparison of the FPRs with some {m,n,a) on the non-negative real axis of g. 

result F{g) very well: the maximal value of the ratio |(p) — F{g))/F{g)\ is 0(10“®). 
Thus we find that our interpolating scheme in this example works quite well at least along 
the positive real axis of Of course this result is not new and has been already seen in the 
previous studies [2, 3]. Here we ask another question. Suppose we perform naive analytic 
continuation of the interpolating function Fm}n{g) to the whole complex plane of g. Then, 
does the interpolating funetion FQ^Q^\g), which is very precise along the positive real g, still 
gives a good approximation beyond the positive real axis? 

In order to answer this question, we plot the quantity® \F^l‘^\g)/F{g) — 1| against \g\ 
for some arg(p) in fig. 3. First, we easily observe in fig. 3 [Left top] that the FPR gives 
very precise approximation for a.vg{g) — tt/IO, whose relative error is at worst. 

Fig. 3 [Right top] shows that this is true also for arg( 5 () = 47r/10, albeit not as accurate as 
that for arg(p) = tt/IO. We therefore conclude that the FPR can give good approximation 
of the exact function even beyond the real positive axis. However, as we further increase 
arg(p) to TT, the approximation starts becoming worse. As seen in fig. 3 [Left bottom], 
the FPR still gives good approximation for a.vg{g) — Gtt/IO but the relative error becomes 
(9(0.1%). Since we have the exponentially suppressed corrections in the weak coupling regime 
for 7r/2 < arg(p) < 37r/4, which comes from the nontrivial saddle points, we recognize that the 
exponentially suppressed corrections are responsible for this error. For 37r/4 < arg(p) < 57r/4, 
the contributions from the nontrivial saddle points become exponentially growing. Since the 
FPR lacks this information, the FPR should show very large error for 37r/4 < a.vg{g) < 57r/4 
in small-I p I regime. Indeed we have (9(100%) error on the negative real axis as seen in fig. 3 
[Right bottom]. In fig. 4, we summarize validity of approximation by the FPR. The shaded 
part shows the region where the FPR has more than 5% relative error. We also draw the 
zeros and poles of the rational function {Fq°^Y associated^ with the FPR, which give branch 
cuts of the FPR. From this figure we observe some points: 

®Note that “ -P(5))/-P(5) does not take real value for complex g in general. 

^ If we consider general FPR f///, then its natural associated rational function is 
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Figure 3: The quantity \{FQ^g‘^\g) — F(g))/F{g)\ is plotted against \g\ for some arg( 5 '). 

• The shaded region looks like a fan, whose radial bounding lines are close to the anti- 
Stokes lines of the small-^f expansion. This is natural since the dominant part of the 
small- 5 ' expansion changes across the anti-Stokes lines and the FPRs do not know this 
information in small-j^l regime. The radius of the fan should be finite since the FPR 
gives the correct large -5 expansion by construction even across the anti-Stokes lines. 

• The boundary of the (shaded) fan-like region is similar to the region surrounded by the 
origin, poles and zeros of the rational function (Fgy^^( 5 ))^. Especially the anti-Stokes 
lines are close to the lines between the origin and the first poles measured from the 
positive real axis®. This would be natural because when the exact function F{g) does 
not have singularities around the poles of (Fg^Q'^^g))'^, then the FPR differs from the 
exact function by a large amount in the neighborhood of its poles. Hence for this case, 
the FPR clearly gives bad approximation around the poles. 

• There is a branch cut from the pole at the negative real axis to 5 = — 00 . Since the 
partition function F[g) has the branch cut on ( 0 , —cxo), we interpret that the branch 
cut of the FPR approximates the “true” branch cut of the exact result. 

• Although the small -5 expansion has the Stokes line at the imaginary axis, the FPR does 
not detect this Stokes line. This is because the dominant part of the small -5 expansion 
does not change across this Stokes line and only the sub-dominant part changes. 

®The first poles are located at arg( 5 ) ~ ±2.358 while the anti-Stokes lines are oriented along arg( 5 ) = 
37r/4~ 2.356. 
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Figure 4: The region where the interpolating function gives bad approximation. In 

the shaded region, the ratio — 1| is larger than 0.05. We also plot “zeros” (the 

symbol “+”) and “poles” (the symbol “*”) of the rational function (Tgy^^(fi'))^ associated 
with the FPR. The blue solid lines denote the branch cnts of 

These results lead us to the conjecture about the general feature of FPR that each branch 
cut of the FPR has the following possible interpretations. 

1. The branch cut is an artifact of the approximation by the FPR. 

2. The physical quantity F{g) has a branch cut near the branch cut of the FPR. 

3. Near the branch cnt, one of pertnrbation series of F[g) changes its dominant part. This 
case implies the following two possibilities on Stokes phenomena. 

(a) We have an anti-Stokes line of the perturbative expansion near the branch cut. 
This possibility occurs most likely for the hrst branch cnt measnred from the 
specific axis where the interpolating function is constructed. 

(b) We have a Stokes line near the branch cut, whose diagonal multiplier is different 
from 1. 

One of immediate questions here is if these featnres are particular for this problem or true 
also for other problems. We will explicitly check this in sec. 4 that this is true also for the 
partition functions of the Od Sine-Gordon model. This seems to hold also in other examples 
such as BPS Wilson loop in 4d W = 4 Super Yang-Mills theory, energy spectrum in Id 
anharmonic oscillator etc [24]. Another important question is if we can construct another 
interpolating function, which gives good approximation in the region, where the interpolating 
fnnction along the positive real axis becomes bad. One natnral way to do this is to construct 
interpolating functions along a specific axis with |arg( 5 ')| > 37r/4, where dominant part of the 
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Im(g) 



Figure 5: Similar plot as fig. 4 for the FPP 

small-p expansion comes from the non-trivial saddle points, and then to extend these to the 
complex g. In next subsection, we will perform this by considering interpolating functions 
along the negative real axis. 

Remarks 

One might ask if the FPRs Fmji with the same (m, n) but different a gave similar results. 
In this example we find that the results strongly depend on a. For instance, let us look at 
fig. 5, which is similar to the plot in fig. 4 but for Note that the FPR for this case 

becomes the FPP [2]. On the real positive axis, this interpolating function gives maximum 
error of about 0.6%. It is easy to see from fig. 5 that the analytic structure of FQ^Q'^^\g) is 
very different from the one of Fg^Q‘^\g). Especially, when we start from the real positive axis 
and go towards the anti-Stokes line at arg( 5 () = 37r/4, we encounter many branch cuts for this 
case. Another important difference is that the first branch cuts from the real positive axis 
are located at the pretty smaller angle arg(p) ~ ±1.21, which is quite far from the location 
of the anti-Stokes likes. This seems to be the reason why the FPP FQ^g‘^^\g) gives the worse 
approximation than Fg^Q‘^\g). 

3.2 Interpolation along negative real axis 

In order to construct another interpolating function precise for |arg( 5 ') | > 37r/4, let us consider 
interpolating function along the negative real axis of g and then consider its naive analytic 
continuation to the whole complex plane. The function F{g) has the following small-^f and 
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Figure 6: [Left] The function F{g) = ie F{g) (black solid), its sniall-(yf expansions 


'<Ns) 


{g) (red dotted) and large-^r expansions Fi^''\g) (blue dashed) for arg( 5 f) = 


TT 


e. 


[Right] The quantity {F^l{t)/F(t) - 1) = (-^l+,TO,n(“^)/-^(“^) “ 1) plotted to t = l^-]. 


large-^f expansions on e-neighbor hood of the negative real axis: 

F{—t + ie) = —s/ieie'^ (2 + 6t^ + 0{t^)) + (2 — 6t^ + 0(t^)) + 0(e) 

= (ir(i/4) - lr(-i/4)r' + ^r(i/4)r^ + o(f-=>)) + o{e), 

F{—t — ie) = +\/:riei^ (2 + + (P(C)) + (2 — 64^+ (!?(('*))+(P(e), 

= (^r(i/4) - lr(-i/4)r' + Lr(i/4)i-2 + o(r=>)) + o(e), 


(3.9) 


with t E IR+. Namely, the dominant parts of the small-expansion and the large-p expansion 
change their signs across the negative real axis. This reflects that the exact function F(g) 
has the sqnare root branch cut on the negative real axis. Instead of F(g), let us consider 
interpolating fnnctions of the quantity: 


F(t) = lim ie '^F(q + ie) 


t G 


g^-t 


(3,10) 


The function F{t) has the small-t and large-t expansions, 

F{t) = 2^/7^ + + 3465\/^i^ + -^ q 


= 


TT 


+ 


Tit 


-1 


Tit 


-2 


Tit 


-3 


V^r (I) 8x/2r (I) 32 x/ 2 r (|) 256 x/ 2 r (I) 


+ 0{t-^) ,(3.11) 


where we have dropped the exponentially suppressed correction 0 {e~'i^) coming from the 
trivial saddle point in the small-t expansion. Denoting interpolating function of F(t) as 
Fm)i(t), we can approximate the original fnnction by using the function 

FiS„j9) = (3.12) 
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Figure 7: [Left] The region where the ratio |-F^+,io,io/-^ “ ^1 larger than 0.05. We also plot 
zeros (+) and poles (*) of the rational function The blue solid lines denote the 

branch cuts of T^+ioio- [Right] Similar plot as the left for 

The interpolating function reproduces the small-p and large-p expansions of F(g) 

for g E M_ ± ie up to certain orders. Indeed F^^^^(g) approximates F{g) quite well along 
the negative real (+ie) axis® as seen in fig. 6 [Right]. Especially lo(fi') has about relative 

0.8% error at worst. 

Let us consider general complex g in the interpolating function Fj^^^(g) as in the last 
subsection. Then, unless we cross the anti-Stokes line or branch cuts, gives the 

correct large-^- expansion and dominant part of small-p expansion. If we cross the anti- 
Stokes line, then the small-^r expansion is dominated by the contribution from the trivial 
saddle point and F^^mnig) should fail to approximate F{g) in small-|(yf| regime. Also, if we 
cross the branch cut particular to the FPR, then the FPR will pick up an extra phase and 
also break the approximation. This extra phase could be trivial depending on the number of 
times, where we cross branch cuts. 

The validity of the approximation by [g) is snmmarized with its analytic property 

in fig. 7. One of important differences from the interpolation along the real positive axis is 
that the (semi-)circular branch cnt of Fl±ioio(fi') snrrounds all the poles and zeros on the 
right plane. Therefore if we go across the circular branch cut on the right half plane, then 
-Tl±ioio( 5') undergoes a sign flip and hence it fails to approximate F{g) across the circular 
branch cut. However, these FPRs also have a branch cut along the positive real axis and 
crossing this branch cut on the positive real axis leads to another flip in the sign. As a result 
^L±io 10 (s') recovers the correct sign and gives the good approximation to F[g) again. That 
is why we see the disconnected nnshaded regions in fig. 7. 

Another important difference is that the first branch cuts as measured from the negative 

® Similar result holds also for 
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Im(g) 



Re(g) 


Figure 8: [Left] First poles and zeros of {Fm,n^y from the negative real axis for various 
{m,n). The symbol “+” denotes zeros and denotes poles. [Right] The region where the 
patch F'{g) of the best FPRs gives more than relative 5% error. 


real axis slightly deviate from the anti-Stokes lines. While the anti-stokes lines are oriented 
along a.vg{g) = ±37r/4 ~ ±2.356, these hrst cnts are located at the angle^° SiYg{g) — ±2.504. 
However, we expect that as increasing the values of {m,n) in the FPR Fm,n \ the first poles 
of {Fmln^Y will approach to the first zeros from the negative real axis and the first cuts will 
finally vanish for sufficiently large (m, n). Indeed we can easily observe in fig. 8 [Left] that 
the pair of first zeros and poles seem to converge to the same point as increasing (m, n). Thus 
we conclude that the first pair of the branch cuts are the artifact of the approximation by the 
FPR with insnfficiently large (m,n). We call this type of singnlarities “fake singularities”. 
The above result would be natural because the FPR with larger (m, n) tends to give better 
approximation in this problem and may improve the validity of the approximation near the 
anti-Stokes lines^^. In next section we will see that FPRs in the Od Sine-Gordon model have 
similar features. 

Finally let us find good approximation of F{g) in region as wide as possible by patching 
the best interpolating functions along the positive and negative real axis. In fig, 8 [Right], 

^*^Note that when we write values of argfg), we always denote the values measured with respect to the 
positive real axis with counterclockwise. 

When either of the small-g or large -(7 expansion is convergent as in this problem, we expect this tendency 
because the convergent expansion itself gives very precise approximation inside its radius of convergence and 
we can regard their FPRs with large (m, n) as analytic continuation of the convergent expansion to whole 
range of values of g. However, it is nontrivial in general whether the decreasing behavior of errors in the 
FPRs will be monotonic or not, although the error of the particular FPR Fm!n^ in this problem seems to be 
monotonically decreasing. For example, in 2d Ising model with finite volume, both low and high temperature 
expansions of specific heat are convergent For this case, error of its FPR with fixed a tends to decrease but 
not monotonically decrease as increasing {m,n) [3]. 
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Figure 9: Summary of Stokes phenomena for the small-^r expansion in the partition f un ction 
of the Od Sine-Gordon model. The blue solid lines denote the Stokes lines while the green 
dashed lines denote the anti-Stokes lines. 


we draw range of validity of approximation by 


ns) 


\ for 

ldS,io(s) for 
id'dbo(s) for 


|arg( 5 f)| < 37r/4 
37r/4 < arg( 5 ') < tt 
—TT < arg(p) < —37r/4 


(3.13) 


This indicates that the patching F'{g) has 5% or better accuracy in the very wide region. 


4 Partition function of zero-dimensional Sine-Gordon 
model 


Let us consider the partition function of the zero-dimensional Sine-Gordon model: 

1 1 • 2 

F{g) = — (4.1) 

\/9 J-tt/2 

which was considered by Gherman-Koroteev-Unsal in the context of resurgence [20]. As in 
the last section, this integral can be evaluated exactly as 

(i) ■ 

The function F(g) has the following small-^- and large-^- expansions 




( y^(8 + 4g + 9g^ + 0(g^)) 

\ \/^(8 + 4g + 9g^ + 0(g^)) - (8 - 4g + 9g^ + 0(g^)) 

{ Vi(^ + 4g + 9g^ + 0{g^)) + {8 - 4g + 9g^ + 0{g^)) 


for 8 ixg{g) = 0 
for arg(fif) G (0, tt) , 

for arg(fif) G {-tt, 0) 
(4.3) 
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Figure 10: [Left] The partition function of the Od Sine-Gordon model (black solid), its 
small-^f expansions Fs^'‘\g) (red dotted) and large-^r expansions (blue dashed) for 

a.vg{g) = 0. [Right] The function {Fm}n/F — 1) is plotted to g for some {m,n,a). 

while large-^f expansion is given by 

F{g) = - ^g-^ + ^g-^ - + 0{g-^)^ . (4.4) 

Notice that the small-(yf expansion has a Stokes line at arg( 5 () = 0, tt and anti-Stokes lines at 
arg( 5 f) = ±7r/2 as summarized in fig. 9. We can again understand this from the viewpoint 
of saddle points analysis. Saddle points of the integration are given by a:* = 0, ± 7 r/ 2 . At the 
saddle points, the action S{x) = ^ sin^ x takes the values 

S{x^ = 0) — 0, S{x^, = ±7r/2) = — 

We can pick up all the saddle points through steepest descent except^^ for arg( 5 f) = 0. 
We again have relative minus sign in the contributions from the non-trivial saddle points 
x^, = ±7r/2 because directions of the steepest descent through x^ = ±7r/2 are opposite 
between the cases for 0 < arg( 5 ') < tt and — tt < arg(p) < 0. As in the Od theory, below 
we consider interpolating functions along the positive and negative real axis, and study their 
analytic properties as complex functions. 


(4.5) 


4.1 Interpolation along positive real axis 


We start with interpolating functions along the non-negative real axis of g and then analyt¬ 
ically continue to complex coupling. For this case, we have the following small-and large- 


expansions 


Ns 


Ff^\g) = ^Skg'", Sk 

k=0 


^ 2'=U(fc + 1/2) 

r(i + i)r2(i/2)’ 


Note that for arg(^) = 0±, direction of the steepest descent around x* = ±7r/2 is arg(x) = T(7r/2 — 0±) 
while the one for = 0 is arg(^) = 0±. 
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* 
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^ Re(g) 


Figure 11: [Left] In the shaded region, the ratio — 1| is larger than 0.05. We 

also plot zeros (+) and poles (*) of The blue solid line denotes the branch cuts. 

[Right] First poles and zeros of from the positive real axis for various (m,n). 

FFis) = 9-^1^ E k<,-\ k = (4^6) 

which are compared with the exact result in fig. 10 [Left], In terms of these expansions, we 
can construct interpolating function Fm}i{g) (see app. C.2 for explicit forms). 

In fig. 10 [Right], we compare some FPRs with the exact result for arg(p) = 0. As in the 
Od theory, we find that these interpolating functions well approximate the exact result 
for positive real g. For example, the best approximation Fg'l^^‘^\g) among the FPRs has 
0(0.1%) error at worst. Next let us consider complex g in the best FPR Fg'l^^^\g). We 
summarize the validity of the approximation with analytic structures of Fg'l^^‘^\g) in fig. 11 
[Left]. The shaded region again starts with the first branch cuts seen from the positive 
real axis, which are given by the first poles and zeros of the rational function (Fg . 

The first poles are located at a-'rg{g) ~ ±0.9348, which is not near from the anti-Stokes 
lines at arg(p) = 7r/2 ~ ±1.571. We again expect that the first branch cuts will vanish as 
increasing (m,n) in Fm,n^‘^\g)- Indeed fig. 11 [Right] implies that the first poles will collapse 
to the first zeros for large (m, n) as in sec. 3.2. Hence we conclude that the branch cuts 
are the fake singularities and the FPR Fm,n^‘^\g) with sufficiently large {m,n) would give 
good approximation in the right-half plane of g. In next subsection, we aim to construct 
interpolating functions approximating the exact result on the left-half plane by considering 
interpolating problem along the negative real axis. 
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Figure 12: [Left] The function F(g) = ie^F{g) (black solid), its small-p expansions Fs^‘’\g) 
(red dotted) and large-^r expansions Fi^^\g) (blue dashed) for arg( 5 f) = tt — e. [Right] The 
relative error (Fmjt/F — 1) is plotted to |p| for real negative g. 


4.2 Interpolation along negative real axis 

Let ns consider interpolating problem along the negative real axis. The function F[g) has 
the following expansions 

F{-t + ie) = ^^(8 - At+ 9t^ + 0{t^)) - + 4t - 91“^ + 0{t^)) + 0{e) 

= + 0 {t-^)^ + 0 {e), 

F{—t — ie) — J '^{8 — At + 9t^ + 0{t^)) + i\l^^e 2 t (^8 + At — 9t^ + 0{t^)) + 0{e), 




(4.7) 


with t G M+. The dominant parts of the expansions change the signs across the negative real 
axis as in the (p'^ theory since F(g) has the branch cut on the real negative axis. Hence it is 
more appropriate to consider the function 


F(t) — lim ie^sF(q + ie) 

' €- 5-+0 ^ ^ 


t G 


g^—t 


(4.8) 


The function F(t) has the expansions, 


F(t) - ++ 


(4.9) 


where we have dropped the exponentially suppressed correction 0{e~^) in the small-p ex¬ 
pansion (see fig. 12 [Left] for comparison of these expansions with the exact resnlt of F{t)). 
Then we can construct the FPR Fmjiit) to interpolate these expansions and approximate the 
original function F{g) by 

d±’,™.„(«) = =F»e"*r,S(-«)^ (4.10) 
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Re(g) 



Re(g) 


Figure 13: [Left] The region where the ratio /F —1\ is larger than 0.05. We also plot 

zeros (+) and poles (*) of the rational function (T^opo*)^- The blue solid lines denote the 
branch cuts of [Right] Similar plot as the left for T^Z^gg^ 

The function reproduces the small-p and large-p expansions of F(g) for g E ]R_±ze 

up to certain orders. 

In fig. 12 [Right] we check that the interpolating functions Fm)n{g) well approximate F{g) 
on the negative real axis. We again summarize the validity of approximation by F^^^ (g) 
and its analytic property in fig. 13. Starting with the negative real axis counterclockwise 
(clockwise), the approximation F^g'^g {g){F^^^g [g)) gets worse across the first branch cut as 
in last subsection, which is located around the line Qxg{g) ~ (—)2.2578 and deviates from the 
anti-Stokes lines. However, according to fig. 14 [Left], the first branch cuts seems to shrink as 
increasing (m, n) in . Thus we expect that the branch cut is the fake singularity 

and the FPR F^!^\{g) {F^y^\{g)) with sufficiently large {m,n) well approximates F{g) 
in the left-top-quarter (left-bottom-quarter) plane. 

Let us patch the best interpolating functions along the positive and negative real axis as 
in the theory. In fig, 14 [Right], we draw range of validity of approximation by 


1 


for 

|arg(5)| < 7r/2 


F'ig) = 

FltS(9) 

for 

7r/2 < arg(5() < tt 

(4.11) 

1 


for 

—TT < arg(5() < —tt/2 



This indicates that the patching F'{g) has 5% or better accuracy in the fairly wide region. 


5 Conclusion and Discussions 

We have studied analytic structures of some interpolating functions and discussed their physi¬ 
cal implications. We have proposed that the analytic structures of the interpolating functions 
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Im(g) 



Figure 14: [Left] First poles and zeros of from the positive real axis for various 

(m, n). [Right] The region where the patch F'(g) of the best FPRs gives more than 5% error. 

provide information on analytic property and Stokes phenomena of the physical quantity, 
which we approximate by the interpolating functions. More concretely, we have mainly con¬ 
sidered the roles of the first branch cuts measured from a specihc axis where we construct 
the interpolating functions. If the first branch cuts are not “fake singularities”, then we 
expect that the cuts approximate those of the exact result or indicate locations of Stokes 
or anti-Stokes lines. When the cuts are fake singularities, we should consider next cuts in 
order to get physical implications. We have explicitly checked our proposal in the partition 
functions of the Od theory and the Sine-Gordon model. This seems to hold also in other 
examples such as BPS Wilson loop in 4d W = 4 Super Yang-Mills theory, energy spectrum 
in Id anharmonic oscillator etc [24]. 

We expect that our result is applicable in more practical problems, where we do not know 
exact solutions. One possible application of our results is to use them to find Stokes behavior 
of the physical quantity by studying analytic structures of interpolating functions. 

We have compared the result of the FPR with the recent result [20] of the resurgence in 
the Od Sine-Gordon model. We have seen that the finite order approximation of the resur¬ 
gence give better approximations than the FPR in the region where the FPR gives relatively 
imprecise approximation while the FPR is more precise in sufficiently strong coupling region. 
This implies that the FPR and resurgence play complementary roles. It will be interesting 
to compare them in more detail. 

One of subtle points in our approach is concerned with the fake singularities. It is unclear 
how we should extend the definition of the fake singularities to more general problem. We 
have defined the fake singularity in our two examples as the branch cut shrinking for large 
(m, n) with fixed a. Since the relative errors seem to monotonically decrease as increasing 
(m, n) in our examples, it would be reasonable to expect that the FPRs with sufficiently large 
(m, n) have larger range of validity in the complex ^r-plane. However, in general problem. 
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FPRs with larger (m, n) do not necessarily give better approximations as seen in [3]. Natural 
extension of the dehnition would be to consider a family of FPRs {FmJ'Jis} satisfying 





^ iris,Us -1 

< 

-F ms-\-i ,Tis-\-i -j 

Fig) 


Ha) 


and dehne fake singularity as a cut, which shrinks as we increase s. It would be interesting 
to pursue this direction further. 

In this paper we have focused on analytic properties of the best interpolating functions, 
which provide the best approximation among the given interpolating functions along a specific 
axis. However, when we do not know the exact results, it is nontrivial which interpolating 
functions give relatively better approximation as discussed in [3]. Although the work [3] has 
proposed a criterion to choose the best interpolating function in terms of two perturbative 
expansions, we need information on large order behavior of the expansions to use the criterion. 
It is nice if we use analytic properties of the interpolating functions to determine the relatively 
better interpolating functions without knowing exact results. In our examples, the exact 
results have the square-root type branch cuts. Interestingly the best interpolating functions 
are at a = ±1/2 and hence also have the square-root type branch cuts. Thus the analytic 
properties of the interpolating functions would be helpful to improve such criterion. 

We have seen that each FPR considered here has its own angular wedge of validity. By 
patching the best FPRs along the positive and negative real axis, we have obtained the 
approximation with larger range of validity than each FPR. It would be interesting if one 
can construct single interpolating function, which gives small-^f and large-^f expansions for 
all SiVg{g). One might think that this was conceptually similar to finding connection formula 
between different Stokes domains in exact JWKB method. However there are some important 
differences. One of such differences is that FPR often still gives good approximation even 
across Stokes line while approximation by FPR necessarily breaks down across the anti-Stokes 
line. For example, in the Od theory, the best FPR along the positive real axis gives precise 
approximation even across the Stokes line at arg(p) = 7r/2 unless we approach to the anti- 
Stokes line at 8ixg{g) = 37r/4. Thus the connection formula in JWKB does not seem to give 
useful hints. Let us see one of main difficulties to construct the single interpolating function 
valid for all arg(p) in the toy example: 

oo oo 

Fig) = Foig) + e~9Fiig), with Foig) = Fiig) = 

fc=0 k=0 

where the summations Fo{g) and Fiig) are convergent, and the exact result of Fig) is given 
by their analytic continuations. Let us consider FPR to interpolate small-(yf and large-(yf ex¬ 
pansions of Fig). If we construct FPR along the positive real axis, then the FPR interpolates 
Foig) and the large-p expansion of Fig). Some part of the information about the lavge-g 
expansion is encoded in Fiig), but the FPR does not have this information in small-^f regime. 
This missing information gives exponentially suppressed error on the right plane of g and 
exponentially growing error on the left plane for small j^-l. Similarly if we construct FPR 
along the negative real axis, then that FPR does not have the information that a part of the 
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large-5^ expansion comes from Fo{g). Of conrse we can find much better approximation in 
this example by separately constructing FPRs to interpolate Fo,i(fi') and large-^r expansion of 
analytic continuation of Fo^i{g). However this information is almost equivalent to have the 
exact result and there is no motivation to perform this procedure. 
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A Comparison with resurgence result 


As mentioned in the main text, the authors in [20] have performed resurgence analysis in 
the Od Sine-Gordon model. It is interesting to compare their result with our interpolating 
functions. Their result is 


F(g) 


‘5arg(c/)‘ho(fi') is ^®‘5arg(g)^l(5') 
‘^arg(p)'ho(5^) T ie ‘Saig(g)^l{g') 


for arg(fif) G (0, vr) 
for argO) G (-7r,0) ’ 


(A.l) 


where 

1 /*+ooe*^ ^ 

50d>o,iO) = - / dt e-?Hd>o,i(t). (A.2) 

9 Jo 

The function Hd>o(t) (Hd>i(t)) denotes analytic continuation of Borel transformation of the 
small-expansion coming from the saddle point a:* = 0 (a:* = ± 7 r/ 2 ), which is explicitly 
given by 

B$o(t) = v^2Ai Q,I;l;2t]), 2 U Q, i; 1;-2t]) . (A.3) 


Although we know all the coefficients of the small-t expansions, let us consider finite order 
approximation of the resurgence result as in [20] and compare this with the FPRs. Namely, 
instead of using B$o,i(t), we terminate its small-t expansion up to and use its 

one-point Pade approximant^^: 


pi2N+l)^t^ = 


i + Ef=i4t^’ 


(A.4) 


in the integration (A.2), which reproduces the small-t expansion up to This pro¬ 

cedure is often called (one-point) Borel-Pade approximation. 

In fig. 15, we compare the result of the resurgence with the FPR Fg for^^ 

By resurgence”, we mean the result obtained by the replacement F$o,i(^) ^ 

Strictly speaking, this is so-called diagonal Fade approximant. 

This choice of arg(^) is due to a technical problem in the integration (A.2). As decreasing arg(^), it 
becomes harder to obtain precise values of the integration. Here we expect that this result for arg(^) = ^ 
is almost the same as for arg(^) ^ 0. 
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Figure 15: [Left] — 1| for arg( 5 ') = ^ [Right] The same plot as the left in a 

different scale. 

in (A.2) and (A.l). We first find that the resurgence result gives very precise approximation 
in weak coupling regime. Furthermore all the results of the resurgence give better approx¬ 
imations in the region where the FPR gives relatively imprecise approximation. As we go 
to the large-1 regime, the approximation by the resurgence becomes monotonically worse 
but the 0(g^^) resurgence still gives reasonable approximation at \g\ = 2. However, if we 
further go to very large-| 5 f| regime, the approximation will breakdown^^. For example, the 
0{g^^) resurgence has about 15% error at \g\ = 100 and 140% error at \g\ = 1000. Hence in 
sufficiently strong coupling regime, the FPR always gives the better approximation than the 
finite order approximation of the resurgence. This result implies that the FPR and resur¬ 
gence play complementary roles in describing the exact function at least in this example. It 
will be interesting to compare them in more detail and other examples. 

B Comments on interpolation in Borel plane 

When small- 5 ' (l^'^ge-^) expansion of the function F(g) is convergent, we expect that F{g) is 
very precisely approximated by the FPRs Fm}i with large m{n). How about the case where 
small -5 expansion is asymptotic but Borel summable? For this case, Borel transformation of 
the small -5 expansion is convergent. If we construct FPR-type interpolating function in the 
Borel plane, then one might expect that the “Borel-FPR” approximates the exact result of 
F{g) very well. However, in this appendix, we discuss that the Borel-FPR gives slightly worse 
approximation than the usual FPR at least for the partition function of the Od Sine-Gordon 
model. We have not understood any clear reasons behind this observation. 

Of course, if we include arbitrarily higher order terms, then the resurgence can provide arbitrarily precise 
approximation. This is also the main difference between resurgence method and our interpolating function 
approach. 
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B.l Method 


Let us construct interpolating function for F{g) in Borel plane. For this purpose, it is 
convenient to introduce following quantities: 

Ns Ni 

F'/-\g) = = 9‘Y1 ^Fig) = = / E ‘“N"’ <B-1) 

k=0 k=0 


where 

d = “1“ n, 6 — FI h. 

The parameter M is a real number satisfying 

a > 1, h>Ns + l. 

The Borel transformation of the small-p expansion is 


Ns 




^g+fc-l 


a r(a + *) 


We also define action of B to the large-g- expansion as 


Ni 


= Y 


- r(6 - k) 




(B.2) 

(B.3) 


(B.4) 


(B.5) 


Then we can construct the FPR-type interpolating function BFm}i{t) in Borel plane, which 
interpolate BFs^"\t) and BFi^''\t) up to and respectively: 

^F^lit) - BF^^^\t) = C»(t“+^*), BF^l{t) - BFl^^\t) = (B.6) 


This implies that the original function F{g) is approximated by 


pBia) 

m^n 



dt e 



(B.7) 


which we call “Borel FPR”. Indeed one can show that the function gives Fs^^\g) in 

small-^- regime and Fi^''\g) in large-p regime up to some orders^®. One of subtle points in 
this approach is that different values of M give different Fm^n^ even if we consider the same 
(m, n, a). In this paper, we will not study M-dependence and fix as M = 10. 


B.2 Comparison with usual interpolating function in Od Sine-Gordon 
model 

Let us compare the Borel FPR with the usual FPR in the Od Sine-Gordon model for di.rg{g) = 

0. In fig. 16 we plot 

7 ^( 20 ) 7 ^( 20 ) — j^B{ —1/2) 

-C s r r r r m^n -C r m,n -C 

F ’ F ’ F ’ F ’ 

If M did not satisfy the condition (B.3), then we could not guarantee this property. 
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Figure 16: Comparison of the Borel FPR Fm^n^ and usual FPR Fmji for {m,n) = 
(1,1), (2, 2), (3, 3), (4,4) with (a, M) = (-1/2,10). 

against g for some values of (m, n). We find that the Borel FPRs give slightly worse approx¬ 
imations than the usual FPRs with the same (m, n, a) at least for these four cases. This 
differs from our naive expectation and we have not found any clear reasons for that. We will 
not address this issue further, but it would be interesting to find some good interpretation 
or modification of our construction to get better approximation than the usual FPR. 


C Explicit forms of interpolating functions 

In this appendix, we write down explicit forms for interpolating functions used in the main 
text. 


C.l Partition function of Od theory 


C.1.1 Interpolation along positive real axis of g 

Fi^F\g) = V^(— .) +l) ^ , Fi\^^\g)=y/ 27 Tr(l/ 4 ) ( 8wgr(l/4)+r(l/4)^+27rr(-1/4) -\ 

O’O Vr(l/4)^ J 1^1 yyj V \647r^g^+87rgr{l/4)^Fr{l/4)^F27rr{-l/4)r(l/4) ) 


6.98929^'^ +7.08691^^+1 


1/6 


F(V"fi9)=2-50663,/-37.9117g2 + 10 1532g+l- 

2,2 y 72.4854^^+43.9117^2^10.1532^+1 


f{Y^^(^)=2-50663( 

F^^^^\g)=2.bmQ2,(^2b.bA99g^ +A8.1779g^ +2,2.1A82g^ F^Qg"^ + 1^ 

f(V^Y^)=2.50663,/ 324.019^3+110.261^2+10.0304^+1 ~ 

3,3 y 619.509^4^420.201^3^116 261^2^10 0304^ + 1 

F(V®h9)=2.50663(-,-28 252592+8.09979+1-^-p/® 

3,3 V 197.465^3+256.834^4+145795^3+46.2525^2+8,0997^+1 y 


(C.l) 
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-1/14 


1/10 


F^Y^^^(^)=2-50663(93.3994^'^+220.976^®+239.216^^ + 155.758^^+42^2^;l) 

jr(V2).^2 50663 / 3224 56^4 ^1303.49^3 +238.239^2 +22 8745^ + 1 

4,4 ■ y 6165.22^5^4576.^4^1440.74^3^244.239^2^22.8745^+1 ’ 

j^(l/10) . X 2 50663 [^._ 14 4369.9^+5.07251.9 + 1 _ 

4.4 ■ [368.8697+752.9549^*+708.68995+403. IO694+ 152.17594+44.436992+5.072519+1 

Y^®^(9)=2.50663(1341.4289^+1038.599®+1475.359'^ + 1294.349®+780.7889^+5949^+5492 + 1^ 

F ^^^‘^ hq ) = l S 12 S ! (9+0-124648)(92+0.1339259+0-01439)(92+0.2038149+0.0156882) _ 

5.5 ■ y (92+0.1442229+0.0144064)(92 +0.2802519+0.0278064)(92 +0.3759059+0.0367401) ’ 

^(1/2). gi 2 S I (92+0.1182579+0.0116795)(92+0.1728699+0.0148819)(92+0.2259079+0.0137177) 

6.6 ' y (9+0.179694)(92+0.1180589+0.011219)(92+0.2343569+0.0211961)(92+0.3229149+0.0291839) ’ 

F^Y^®^(9)=2.50663 (^4562.589^^+20047.39^2 ^42029.49^^+55850.99^°+52692.49^+37452.19®+20758.69'^+9156.79® 


+3248.639®+929.3149^+213.0789®+38.491492+5.282929+1^ 


1 


(C.2) 


C. 1.2 Interpolation along negative real axiapf 9/3\2 \ /oO \ 1/2 

2.1 r 1 j.(i/2).,. „ j +"j+"^U) 

0)0 2 3 \ 2 ’ 1)1 / e \ / r / o \ 4 / o \ 2 r\ / o \ 3 

■- ^^r(|)(^64t2r(|) +8,rtr(|) +7r2j+2,rr(|) 


\ 8tr(|)+,r’ 


?<Y^>(t)=6.28319./ 0.0832t2+0.0357t+0.00906 (t) = i.81./ ^ , (t+0.398)(t2+0 0631t+0.0634 ) 

,2 V t6+o.0907t2+o.ll21t+0.0285 3,3 y (t2 _o.240t+0.105) (t2+o.363t+0.0( 


2 ,, 

^(1/2 ).^^^f _ 

4,4 • V t5+o.06t4+o.081t3+o.011t2+o.0040t+0.00116 ’ 


0626) ’ 


t4+0.40t3+0 069t2+0 015t+0.0044 
06t4+o.081^3+0.011^2+o.0040t+0.1 

F^^/‘^ht) = l 81 I (t+0.30)(t2-0 19t+0.044)(t2+0.31t+0.050)~ 
5,5 ^ ^ • y (t2_o.30t+0.063)(t2+o.050t+0.076)(t2+o.34t+0. 


037) ’ 


(’-+49=1.8 


' 6,6 


(t2_0.23t+0.037)(t2+0.10t+0.036)(t2+0.51t+0.07) 
(t+0.18)(t2_o.30t+0.050)(t2_o.042t+0.07)(t2+o.20t+0.042) ’ 


F^^/‘^\t)=2 ! (^+0-3)(t^+0-003t+0.04)(t2+0.4t+0.05) 

7,7 Y (t2_o.lt+0.06)(t2+o.lt+0.05)(t2+o.4t+0.03) ’ 


?(l/2), 
" 8,8 


(t2_o.3t+0.03)(t2_o.04t+0.03)(t2+o.2t+0.02)(t2+o.5t+0.06) 
(t+0.16)(t2_o.3t+0.04)(t2_o.l4t+0.05)(t2+o.09t+0.05)(t2+o.2t+0.02) ’ 


~(l/ 2 ).^._^ o I (t+0.25)(t2_o.25t+0.026)(t2_o.lt+0.028)(t2+o.09t+0.026)(t2+o.4t+0.05) /p 

9^9 ' Y (t2_o.27t+0.03l) (t2_o. 16t+0.04) (t2+o.019t+0.05) (t2+0.13t+0.034) (t2+o.35t+0.033) ’ ^ ‘ 7 


C.2 Partition function of Od Sine-Gordon model 


C. 2.1 


-1/2 


(2+7r)9+7r-l J 


57r+57r2)93 + 187r(9+87r)92+37r(377r-28)9+27r(10-357r + 167r2) ^ 


^2,2 A/27r^ 16(8+57r+57r2)92 + (_64+4l7r+327r2)9+327r2-707r+20 


^(-1/2) g i4i59( -0-03955949+0. 111122)(92+0.8284939+0.228908) 


' 3,3 


(9+0.26163) (92+0.02730419+0.0618949) 


^(-1/2). ,3 14159/^ (9+0.403343)(9^-0.108029+0.0594072)(92+0.4740719+0.13076) \ ^ 

4.4 (92-0.08795239+0.0413742) (92+0.3573469+0.0482101) J 

^(-1/2) 14159/^ (9^-0.1286459+0.0352297)(9^+0.2626479+0.079582)(92+0.6233639+0.110622) \ 

5.5 (9+0.189845)(92-0.1277729+0.0289318) (92+0.1952929+0.0359479) J 


- 1/2 


^(-1/2)^^^^3 141591 (^+0-292344)(9^ -0.1326069+0.0227078)(9^ +0.1395999+0.0510561)(9^ +0.4498869+0.0759508) 


^6,6 


(92-0.1371829+0.0208248) (92+0.08504849+0.0282403) (92+0.3013569+0.0278662) 


- 1/2 


^(-1/2) 11159 f (^^-0-1301399+0.0157197)(9^+0.06681579+0.0339803)(9^+0.3171799+0.0532544)(9^+0.4894999+0.0644541) ^ 

' (9+0.148792) (92-0.1344879+0.0153154) (92+O.OI7OI889+O.O230983) (92+0.2120319+0.0221752) 


^7,7 


-(-1/2) , / (9+0.228889)(92 - 0 . 1247239+0.0115314)(92 +0.02179099+0.0231509)( 9 ^ +0.2216489+0.0382247)( 9 ^ +0.3913529+0.0486489) 

c ( 9 )—3.14159 - ___ — _\ / o —TTrTTTT:::::::— _/ o — _/ o __ — _ - 


' 8,8 


(92-0.1274969+0.0115172) (92-0.02190339+0.0192889) (92+0.1361239+0.0181046) (92+0.2522339+0.0179858) 


- 1/2 


(C.4) 




1,1 


3^/^\t)=3.14159 


(2+'7r)t+'7r —1 


3ng^i;iegatiy^_^gal 




^+57r^)t® + 187r(9+87r)t^+37r(377r-28)t+27r(10-357r + 167r^ ) ^ 


16(8+57r+57r2)t2 + (_64+4l7r+327r2)t+327r2-707r+20 


(t^-0.0395594t+0.111122)(t^ +0.828493t+0.228908) ^ 
(t+0.26163) (t2+o.0273041t+0.0618949) 


^(-1/2) 14159/^ (t+0.403343)(t^ -0.10802t+0.0594072)(t^ +0.474071t+0.13076) \ 

4,4 (t2_o.0879523t+0.0413742) (t2+o.357346t+0.0482101) J 


- 1/2 


(C.6) 
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14159 

14159 

F^“^/^\t)=3.14159 
Fg^g^/^\t)=3.14159 


(t^-0.128645t+0.0352297)(t^ +0.262647t+0.079582)(t^+0.623363t+0.110622) ^ 


(t+0.189845)(t^-0.127772t+0.0289318) (t^+0.195292t+0.0359479) 

(t+0.292344)(t^-0.132606t+0.0227078)(t^+0.139599t+0.0510561)(t^+0.449886t+0.0759508) \ 

(t^-0.137182t+0.0208248) (t‘'^+0.0850484t+0.0282403) (t^+0.301356t+0.0278662) j 

(t^-0.130139t+0.0157197)(t^ +0.0668157t+0.0339803)(t^+0.317179t+0.0532544)(t^+0.489499t+0.0644541) \ 

(t+0.148792)(t^-0.134487t+0.0153154) (t^+0.0170188t+0.0230983) +0.212031t+0.0221752) j 

(t+0.228889)(t^ -0.124723t+0.0115314)(t^ +0.0217909t+0.0231509)+0.221648t+0.0382247)+0.391352t+0.0486 
(t^-0.127496t+0.0115172) -0.0219033t+0.0192889) (t^+0.136123t+0.0181046) (t^+0.252233t+0.0179858) 


- 1/2 


(C.6) 


C.2.3 Borel-FPR 


^^(-1/2) .^, 256^ / 89260141168201575t^ +753534959720857600t+5438758396457123840 ^ 

1,1 ‘ 2 ^ 3 b (2416615932t+10130476945) J 


- 1/2 




BF^ = 


BF^ g^/^\t)=0.0000263352 


tl^(679261120778364683648t^+2276065430590939023520t+7445431147866433926325) 

. - 1/2 

_^ ^ 

tl^ (t+2.86868) (t2+0.141631t+7.48056) 


BF^ 4^/^\t)=0.0000263352 


(t+4.2183) (t^-1.833t+7.96513) +4.69262t+14.8922) \ 

ti«(t2 i.23764t+5.48437)(t2+4.06556t+6.27685) 


- 1/2 




0 ' 
(C.7) 
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